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■ Abstract. A full treatment of decoherence and dephasing effects in BEC 
, interferometry has been developed based on using quantum correlation functions 

■ for treating intcrferomctric effects. The BEC is described via a phase space 
^ I distribution functional of the Wigner type for the condensate modes and the 

^ ■ positive P type for the non-condensate modes. Ito equations for stochastic 

, condensate and non-condensate field functions replace the functional Fokker- 

Cs| ■ Planck equation for the distribution functional and stochastic averages of field 

I function products determine the quantum correlation functions. 

o : 

^ ! 1 Introduction 

' Bose-Einstein condensates (BEC) in cold atomic gases are an example of a 

quantum system on a macroscopic scale. Well below the transition temperature 
essentially all bosons occupy small number of single particle states (or modes) 
- in simple situations only one mode. Interferometry based on BECs (such 
I as by splitting a BEC in a single trap into two traps and then allowing the 

O^' BECs to recombine) offers possible improvements in precision over single atom 

^ I interferometry by a factor given by square root of the boson number N [1, 

2]. Such BEC interferometry is based on having almost all bosons in one (or 
perhaps two) modes. 

^ . A typical double-well interferometry experiment involves starting with a 

" " ' BEC in a single well trap, then changing the trap to a possibly asymmetric 

double-well and then back to a single well. Asymmetry could lead to excitation 
of bosons to higher energy states of final trap, or to changes to spatial interfer- 
ence patterns. The process of exciting one boson from the ground to the first 
excited state for a single well trap via two quantum pathways is shown, both 
pathways involving an intermediate double well trap. The near degeneracy of 
energy levels for the intermediate asymmetric double well facilitates the boson 
transfer to the excited state. The two non-observed quantum pathways involve 
the boson transfer occurring in different halves of the process, and superposition 
of the two quantum transition amplitudes can lead to interference effects in the 
excitation probability. 

Interferometry experiments involving measurements of boson positions are 
best described in terms of quantum correlation functions, which are expectation 
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values of products of bosonic field operators specified at different spatial points, 
and are related to many-boson position measurements [3]. If boson-boson in- 
teractions were absent and the BEC isolated from the environment, idealised 
forms of quantum correlation functions would result, with interferometric effects 
clearly visible. Interactions of the BEC with the external environment (such as 
fluctuating trap fields) and internal boson-boson interactions tend to degrade 
the interference pattern. Boson-boson interactions can result in dephasing (as- 
sociated with interactions within condensate modes) and decoherence effects 
(associated with interactions causing transitions from condensate modes) even 
when external environmental effects arc absent. 

Recently a simple theory of double- well BEC interferometry has been devel- 
oped [4] based on a two mode approximation and allowing for possible fragmen- 
tations of the original BEC into two modes - which may be localised in each well. 
Many previous theories ignore fragmentation, with all bosons occupying a single 
condensate wave function which satisfies the standard Gross-Pitaevskii equation 
(sec [5] and references therein). The two-mode theory was developed from the 
quantum principle of least action and gives self consistent coupled equations 
for the two mode functions (generalised Gross-Pitaevskii equations) and for the 
amplitudes (matrix equations) describing fragmentation of the BEC. Self con- 
sistency results in the mode functions depending on the relative importance 
of various ways the BEC can fragment, whilst the fragmentation amplitudes 
depend on the mode functions. Numerical studies are planned. However, only 
transitions within the two condensate modes are allowed for, and although some 
dephasing processes are included, decoherence processes associated with elemen- 
tary collective excitations (Bogoliubov) and single boson excitations (thermal 
modes) are not. The theory is also restricted to small N. 

A full treatment of decoherence and dephasing effects on quantum correla- 
tion functions is required, and this is the subject of the present paper. Sev- 
eral possible approaches to developing such a theory may be identified. These 
all have desirable features such as exploiting the physics of large occupancy 
differences for condensate and non-condensate modes, not being restricted to 
small boson numbers and avoiding the explicit consideration of large numbers 
of modes. Collective and single boson excitations from condensate modes, and 
possible fragmentation effects can be allowed for, though the presence of these 
processes may not be explicit. 

One such approach is a master equation method [6], in which a conden- 
sate density operator is defined for which a master equation is derived allowing 
for interactions with non-condensate modes, which constitute a reservoir. The 
difficulty with this method is that it is hard to evaluate the non-condensate 
contributions to quantum correlation functions. A second approach is based on 
a Heisenberg equation method that has been applied in numerous many-body 
theory cases. Heisenberg equations for field operators and products of field op- 
erators are derived, and taking the expectation values with the initial density 
operator results in a heirarchy of coupled equations for quantum correlation 
functions. An ansatz (such as assuming that a suitable high order correlation 
function factorises) produces a truncated set of coupled equations from which 
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correlation functions of the required order can be calculated. The problem with 
this method is that it is hard to confirm the validity of the ansatz. 

The present approach is a generalised phase space method, but involving a 
distribution fmictional rather than a distribution function [7]. The field opera- 
tor for the bosonic system is written as a sum of condensate and non-condensate 
mode contributions [8]. The BEC state is described by a density operator which 
satisfies the Liouville-von Neumann equation (LVN) , and which is mapped onto 
a phase space distribution functional. The latter has the feature that the highly 
occupied condensate modes are described via a generalised Wigner representa- 
tion (since the bosons in condensate modes behave like a classical mean field), 
whilst the basically unoccupied non-condensate modes are described via a pos- 
itive P representation (these bosons should exhibit quantum effects). The LVN 
equation is replaced by a functional Fokker-Planck equation (FFFE) for the 
distribution functional, which is based on the truncated Wigner approxima- 
tion [7] that can be applied when large condensate mode occupancy occurs. 
The FFPE are finally replaced by coupled Ito stochastic equations (c-number 
Langevin equations) for condensate and non-condensate field functions, where 
the Ito equations contain deterministic and random noise terms - identifiable 
from the FFPE. Stochastic averages of the field functions then give the quan- 
tum correlation functions. There are no obvious difficulties with this approach, 
though further development will be needed to explicitly incorporate Bogoliubov 
collective excitations. 



2 Generalised phase space functional theory 

The Hamiltonian in terms of bosonic field operators \E'(r), ^'(r)^ is 
H = /dr(|-V$(r)t.V$(r)-F$(r)V$(r)) 



/ 



dr|$(r)^$(r)t$(r)$(r), (1) 



where the boson mass is m, the trap potential is V{Y,t) and g specifies boson- 
boson interactions in the zero range approximation. 

The field operator 'J'(r) is written as sum of a condensate term ^c(r) and 
a non-condensate term ^jvc(r) [6] 

n 

$c(r) = ci(/>i(r) +C2</)2(r) $jvc(r) = ^ c^^), (2) 

which are defined in terms of n orthonormal mode functions (r) and bosonic 
mode annihilation operators Cfc. Condensate and non-condensate contributions 
to the field operators commute. 

Replacing ^'(r) by 4'c(r)-|-\E'Arc(r) the Hamiltonian is sum of three terms [6], 
Hamiltonians He and H^c for the condensate and non-condensate - which are 
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of the same form as in Equation ([T]) , and the interaction V between condensate 
and non-condensate. This is the sum of three contributions, which are hnear, 
quadratic and cubic in the condensate operators. 

Spatial coherence effects for interference experiments in BECs may be de- 
scribed via quantum correlation functions, which depend on the density operator 
p for the bosonic system and are defined by 

G^(ri,r2, ..,rjv;sjv, ..,83,8^ 
= Tr(p(t)$(ri)t..$(rA,)t$(sAr)..$(si)). (3) 

The quantum correlation function with — Si{i ^ l,...,N) determines the 
simultaneous probability of detecting one boson at ri, .., the A^th at tat, see 
[3]. It is evident that the quantum correlation functions will contain conden- 
sate terms (describing the main interference effects) , non-condensate terms and 
mixed terms involving both condensate and non-condensate operators (describ- 
ing effects degrading the interference patterns). 

In the phase space functional method the density operator p is first mapped 
uniquely onto a characteristic functional xKc(i")' Cc('^)' ^Afc('^)' ^ivc('")] °^ ^^"^ 
four functions ^J(r), ^c'(r), 'f^c('") ^^d ^Nci^) 

= rr(pexpz J drUc(r)*^(r) + $c(r)eJ(r)} 

xexpi J dr{^^c(^)^Nci^)} exp« J dr{$wc(r)^^c(r)}) (4) 

The characteristic functional is of the Wigner (W) type for condensate modes 
and the positive P (P+) type for the non-condensate modes. 

The (quasi) distribution functional P['ilj(y{r),tp'^{r), ■(/'^^(r), i/'jvc('")] involves 
four field functions ^(y{r),tjjQ{r),tl;jy^{r),ip~l^fj{r) corresponding to field oper- 
ators 5'c(i"), 5'c(r)^ ^'jvc(r) and \l/jvc(r)^. Although non-unique and possibly 
negative (and hence not interpretable as a probability distribution) , it is required 
to determine the characteristic functional x[Cc('')i 'Cc('')' 'CArc('')i ^JVcC*")] ^ 
functional integration process with weight function w(V'i, , ■■, ^i, ipt ^ ••> V'n ; 
given by J|(Ar,,) 

i 

xKc(r),^J(r),^^c(r),^^c(r)] 

i?Vc D^^i ^Vatc D^i^+c P[^cir), V'J(r), V'jvcW, V'^c W] 

xexp*! drUc(r)^J(r)-FVc(i-)eJ(r)} (5) 

xexpi / dr{f^c'(i')V4c('")}expi / dr {V'A^c(I■)^^c(l■)}• 



Quantum averages of symmetrically ordered products of condensate field 
operators {4'p(ri)....\I'^(rp)4'c'(sg)..\I'c(si)} and normally ordered products of 
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non-condcnsatc field operators '$\^(j{ui)'^]^(j (u2 ) . . . .4' jy^ (u^ ) nc (vs ) ■■'^nc (vi ) 
are given by functional integrals of the distribution functional P[tpc, tpQ, V'tvc] 
with products of field functions,. where the condensate field operator ^'c(ri)^ is 
replaced by ij)Q{ri), ^c{^j) is replaced by ip{sj) and with analogous replace- 
ments for the non-condensate field operators. 

Tr[p{%{v^)...%{v^)^c{^,)..^c{^i)] (6) 

X$]y(^(Ui)...$jvc(Ur)*ArcK)..$jVc(vi)] 

= //// ^Vc D^i'i D^NC D'i^tfc Pii^ci^UU^), VivcW, ^^^c (r)] 

xV'^(ri) ..V'c(rp) tPcisq) -^0(81) V'^c(ui) -V'^cK) Vjvc(vs)-^jvc(vi)- 

Symmetric ordering is defined as the average over all N{p, q) = {p + q)\ permu- 
tations of the factors \l/^(ri)..\E''''(rp)4'(sg)..4'(si). These results plus equal time 
commutation rules give the quantum correlation functions. 

For example, the first order quantum correlation function that is used to 
exhibit macroscopic spatial coherence in BECs is given by 

Gi(ri;si) = /$(ri)t$(si) 



-ls{r,-s,) (7) 

D^^jc D^i4j D^^NC D^^NC ^[V^cW, V'J W, V'ATCW, V'^cW 

x(V'J(ri) + ip+c{ri)){^ci^i) + V'jvc(si))- 

The result includes pure condensate terms, pure non-condensate terms and 
mixed terms. The delta function term arises from the difference between normal 
and symmetric ordering for the condensate terms. 

The Liouville-von Neumann equation for the density operator replaced by 
the functional Fokker-Planck equation for the distribution functional by using 
the correspondence rules involving functional derivatives, such as ^'(^(s)^ ^ 

(ijjcis) + ^S/5^/j+{s)) P[i/jc, ..,yj%c] and ^Ncis)?^ {^P^ci^)) Pi^^c^ -^^nc]- 
The density operator is multiplied by various field operators in the LVN equation 
and the overall effect on the distribution functional is obtained by applying such 
rules in succession. The general form for the FFPE after applying the truncated 
Wigner approximation to remove the third order functional derivatives is 

d 

f ^ 5 

^ a,/3=l ^'^"^^^ 

xP[V'c>^i'C'^Jvc>^^c]' (8) 
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where (j^i = 'ip(j,4>2 = '^ci'^a = '^NC^'t'i = ''P'nc- drift vector Aa involves 
the fields ipQ^ip^^'ipj^Q^'ip'^fj and their spatial derivatives. The positive definite 
diffusion matrix D^p also involves these fields. 

The functional Fokker-Planck equation for the distribution functional is then 
replaced by Ito stochastic field equations for ipf^, -0^, "f/'Arc V'jvc which are now 
regarded as time dependent stochastic fields. The general form for the stochastic 
field equations is 

d 

■jrAa = Aa{ii^c^ V-c, V'atc. ^nc) + XI '^"/3(^C: i'c^ V^wc)r/3(r, t). (9) 

p 

The first term is deterministic and is associated with the drift vector. The 
second term is a random noise term and is associated with diffusion matrix, the 
matrix d being related to the diffusion matrix via D = d<f^ . The r^(r,t) are 
Gaussian-Markov random noise terms 

r„(ri,ii) =0 r„(ri,ti)r^(r2,i2) = <5a^<5(ti -t2)<5(ri -r2) (10) 

where the bar denotes a stochastic average. 

Quantum averages of the field operators are now given by stochastic averages, 
which are equivalent to the functional integrals in Equation ([6]) involving the 
distribution functional. 

= V'c('"l) ■•V'c('"p) V'cK) •■'0c(si) V-^cK) ••V'wcK) V'ArcK)-^Arc(vi) 

(11) 



3 Conclusion 

It is shown how the quantum correlation functions required for describing in- 
terferometry using BECs can obtained via stochastic averages of products of 
field functions, which satisfy Ito equations derived from the functional Fokker- 
Planck equation for the phase space distribution functional that represents the 
quantum density operator. The phase space distribution functional is of the 
Wigner type for the condensate modes and the positive P type for the non- 
condensate modes. Decoherence and dephasing effects in BEG interferometry 
with large boson numbers are fully treated, unlike previous theories in which 
only condensate modes are considered or the theory is restricted to small boson 
numbers. 
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